ANALYSIS OF THE QUASI-NONLOCAL APPROXIMATION OF LINEAR AND 

CIRCULAR CHAINS IN THE PLANE 



PAVEL BELIK AND MITCHELL LUSKIN 

Abstract. We give an analysis of the stability and displacement error for linear and circular atomistic 
chains in the plane when the atomistic energy is approximated by the Cauchy-Born continuum energy and 
by the quasi-nonlocal atomistic-to-continuum coupling energy. We consider atomistic energies that include 
Lennard-Jones type nearest neighbor and next nearest neighbor pair-potential interactions. 

Previous analyses for linear chains have shown that the Cauchy-Born and quasi-nonlocal approximations 
reproduce (up to the order of the lattice spacing) the atomistic lattice stability for perturbations that 
are constrained to the line of the chain. However, we show that the Cauchy-Born and quasi-nonlocal 
approximations give a finite increase for the lattice stability of a linear or circular chain under compression 
when general perturbations in the plane are allowed. We also analyze the increase of the lattice stability 
under compression when pair-potential energies are augmented by bond-angle energies. Our estimates of the 
largest strain for lattice stability (the critical strain) are sharp (exact up to the order of the lattice scale). 

We then use these stability estimates and modeling error estimates for the linearized Cauchy-Born 
and quasi-nonlocal energies to give an optimal order (in the lattice scale) a priori error analysis for the 
approximation of the atomistic strain in (? E due to an external force. 



1. Introduction 

The quasicontinuum (QC) method |23j is a general approach for coupling atomistic models with Cauchy- 
Born continuum models to achieve both accuracy and efficiency. Many authors have improved, extended, 
and analyzed the QC method and related atomistic-to-continuum coupling methods |^l^[5ll MI^[T^[T5l[17r - 

iiipaii)]. 

In this paper, we give a linearized analysis of the stability and strain error in for linear and circular 
atomistic chains in the plane when the atomistic energy is approximated by the Cauchy-Born continuum 
energy and by the quasi-nonlocal atomistic-to-continuum coupling energy |28j . We consider first atomistic 
energies that include only Lcnnard- Jones type nearest neighbor and next nearest neighbor pair-potential 
interactions, and we then consider atomistic energies that also include bond-angle interactions. 

We chose the quasi-nonlocal atomistic-to-continuum coupling energy because uniformly spaced linear 
chains (one-dimensional lattices) are equilibria for the coupling energy just as they are for the atomistic and 
Cauchy-Born continuum energies [Tll2"l] . This property is called patch test consistency. Patch test consistent 
extensions of the quasi-nonlocal energy to finite range interactions have been given in [9l ll6i[27] . 

A uniformly strained one-dimensional lattice modeled by a Lennard-Jones type atomistic interaction loses 
stability when the strain reaches a critical value (the critical strain) . We seek to estimate the critical strain 
for quasicontinuum energies and to then compare them with the critical strain for the atomistic energy. We 
define such estimates to be sharp if they are exact up to the order of the lattice spacing. Sharp lattice 
stability and optimal order (in the lattice scale) strain error analyses of the one-dimensional quasi-nonlocal 
approximation have been given in [6l[7l l21[[24] . In this paper, we give a sharp stability and optimal order 
strain error analysis of the linearized problem for the quasi-nonlocal approximation of some simple objective 
structures that are generated by a single affine mapping |13j . Such objective structures include linear chains, 
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circular chains, and helical chains. Wc focus on linear and circular chains in the plane, and we note that our 
analysis allows general planar perturbations. 

Previous analyses of linear chains have shown that the quasi-nonlocal approximation reproduces the 
lattice stability for perturbations that are constrained to the line of the chain. However, chains can undergo 
buckling instabilities under compression when general planar perturbations arc allowed. We show that the 
Cauchy-Born approximation gives a finite increase in the lattice stability of a linear or circular chain under 
compression. We also analyze the increase of the lattice stability under compression when pair-potential 
interactions are augmented by bond-angle energies. 

We restrict our analysis to the classical quasi-nonlocal approximation of chains with next nearest neighbor 
interactions |28j . The analysis in [16] can likely be utilized with the analysis in this paper to obtain a quasi- 
nonlocal analysis of linear and circular chains for finite range interactions, but the arguments would be 
considerably more technical than the analysis presented in this paper. We give our error analysis for the 
linearization about linear and circular chains. The linear analysis in this paper can also likely be extended 
to a nonlinear analysis by utilizing the inverse function theorem techniques developed in [24) , but the details 
of this analysis would also greatly increase the complexity of the analysis that wc present. For simplicity and 
clarity of exposition, we will present our results and analysis for chains in K 2 , but we note that our analysis 
generalizes directly to helical chains in R 3 . 

In Section [2j we define the energy of a chain of atoms and its Cauchy-Born approximation. In Section G2 
detailed derivations of the first and second variations of the atomistic and Cauchy-Born energies are pre- 
sented. We give sharp lattice stability results for unconstrained periodic chains in Section [4] and for circular 
chains in Section [5] We introduce a model for bond-angle energy in Section [6] and derive results for its 
contribution to lattice stability. 

In Section we give estimates for the modeling error due to the Cauchy-Born approximation, and we 
then give an error analysis of the linearized problems in Sections [5] and |H] based on our stability and modeling 
error analyses. Finally, in Section 1 1 01 we define the quasi-nonlocal approximation for periodic chains and 
give a sharp stability and error analysis for its approximation of the atomistic model. 

Wc summarize the results in Section 111! 

2. Definitions 

We will consider chains of atoms y = {ye}g£z, where yi £ R 2 denotes the position of the £-th atom in a 
plane. We will assume that in the reference configuration the chain is straight and that the distance between 
neighboring atoms is e > 0. We will only focus on periodically repeating chains with period N £ N such that 
Ne = 1, where by periodically repeating we mean that ye+N = ye + (un — Ho) for all I £ Z. This definition 
allows us to treat closed chains of N atoms (for which y^ = y^), and also chains that are not closed, but 
such that the shapes of the overall configurations repeat every N atoms (such as straight chains). In 3-D, 
periodically repeating chains would include, for example, helical chains. 

Thus, we can define the space of iV-pcriodically repeating chains, y, by 

y = {y = {yi}i& ■ »eK 2 and y £+N =yg, + (y N - Vo) for all leZ}. 
We will also use the space of iV-periodic mean-zero displacements, hi C y, defined as 

N 

hi = {u = \ui\i<^i : U£ £ M 2 and ue+N = ue for all £ £ Z, and ui = 0}. 

t=i 

We will refer to such displacements as 2-D displacements. 

Finally, we will also consider the subspaces y C y and U C U of one-dimensional periodic chains and 
displacements, respectively, defined as 

y = {y = {yth& e y : » • (0, 1) = for all I G Z} 

and 

U = {u = {u e } e€ z € hi : Ui ■ (0, 1) = for all £ £ Z} 
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when studying linear chains of atoms constrained so that the atoms can only move in the direction of the 
chain. 

For y G y, we define the backward finite differences 



, ye ~ ye-i 
Vi = 



<J< 



i't-v'i-x y't-iy't-x + v'i- 



(n-l) (n-1) 

(«) _ Vi - Vi-i 



for n > 2, 



e ee* 
and write y' = {y' e }e £ z, y" = {y"}e £ z, etc. 

For v, w € IR 2 , we will write v ■ w for the usual inner product in M 2 and \\v\ 
norm. For v,w £ y, we define the inner product 

TV 

(v,w) = e^2ve ■ W£, 



Iv.v for the usual Euclidean 



and the associated norm 



1=1 

f TV \V2 

e vt ■ v i 1 

V 1=1 / 



/ TV 

\\v e \ 

\ 1=1 



1/2 



When estimating approximation errors, we will use the negative- order norm 

(v,w) 



v * 



sup ., ,,, . 

oeu\{o} \\w'\\e e 

For a chain y € y, we can now define the "atomistic" energy, £ a (y), via nearest neighbor and next nearest 
neighbor interactions 

TV 



hi - ye-i\\ 



he+i - ye-i\\ 



TV 



y>i\\vt\\) + 4>(h'i+i + vi\\)], 



(2.1) 



i=i 



and its local Cauchy-Born approximation p][6j[7l|3T] , £ (y) 



N 



£ c *{y) = eJ2£? B (y), 



i=i 



where 



£? B (y) = ^ 



ht+i - yd 



hi - ye-i\\ 



ht+i - yd 



-j> 2 



\ye - yt-i\\ 



In the above approximation, we took into account only the nearest neighbor on either side of an atom at 
the position yi and extrapolated linearly to approximate the position of the next nearest neighbor on either 
side. Thus, the distance between this approximating next nearest neighbor and the original atom is twice the 
distance between the nearest neighbor and the original atom. Note that due to the chain being periodically 
repeating, we can write 

TV 

£ CB (y) = eJ2^(h'i\\) + ^m\y'M- (2-2) 

t=l 

3. Variations 

In the following sections, we will study the linearizations of the atomistic energy (|2.ip and Cauchy-Born 
energy (|2.2[) about periodically repeating deformations y € y satisfying the particular constraints outlined 
below. We begin by deriving expressions for their first and second variations. We first define the strains 

_ M „.,|| and c, ht+i- yt-i\\ _ h'i+i +y'A\ 



Fi = 



\y'e\ 



F 2 = 



\yi+i - yt-i\ 

2e 



(3.1) 
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and we will assume that both F± and F 2 are independent of i. This assumption allows us to study, among 
others, linear (1-D) and circular deformations y £ y with uniformly spaced atoms. Note that a 3-D helical 
chain also satisfies this assumption. 

The first variation of the atomistic energy (|2.1[) about a periodic configuration y £ y satisfying (|3.1|) for 
all I £ Z is 



5£^y)[u}=eJ2 
1=1 

N 

E 

N 

E 



= e 



^y^ + ^ + '? U + i + y^M + i + 4) 



\\Vt+i + Ve\\ 



Fi 



2F 2 



= e 



U'(F 1 ) | 2<P'(2F 2 ) \ 
\ F\ F 2 ) 



I I 2 

Vi-u t -e 



2F 2 



(3.2) 



where we have used the identity (with A = I 2 , the 2x2 identity matrix) 

(y'e+i +V'i)-A (u' e+1 + u' t )=2 y' l+1 ■ Au' l+1 + 2y' r Au' e - e 2 y'l +1 ■ Au'[ +l 

and rearranged the sum using the periodicity of the chain. 
The second variation of the atomistic energy (|2.ip is 

i ( ,„ nl , , ,, h fiM±i + ygy \ (u' e+1 + u' t ) ■ (y' t+1 + y[) (y' e+1 + y[) ■ (v' e+1 + v' t ) 

+ V + M + i + y't\\ J H^ZW KTHfli 
+ W±i±jP K+1 + ^.^ +1 + ^ 



N r 



E 



\\y't+i + y'i\\ 



(3.3) 



0"(2F 2 ) J P, + ^^(/ 2 -P,) 



^r^'iF^ + ^-ih-P^v', 

+ «+i + «J) • U"{2F l )P l + ^§^{h Pi) 



(3.4) 



+ K+i +«/)■( (0'W) - <t>"{2F 1 ))P t + 



{^2F 2 2FT) 



where we have used the projection operators 



Pl = jL®jL and P t - y '^ +y ' e - y '^ +y * 



WvtW WvtW W t+ x + y'A\ Wvi+i + vl\\ ' 

It is easy to see that Pf = Pt and Pf = Pt for all I. If we now apply identity (|3.2[) to the term (u' e+1 + u' e ) 
U"(2Fi)PA (v' t+1 + v[) in (E3), we obtain 



5 2 £*(y)[u,v] = eJ2 U ■ U"{F 1 )P i + Pi)) v't 

t = i l \ 1 / 
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+ K+i + u' t ) ■ {^§^{h - ft)) (v' e+ i + vl) - e 2 u'l +l ■ (V'(2ft)ft) v'l +1 
+ 2u' e+1 ■ (>(2ft)ft) v' i+1 + 2u' r (>(2ft)ft) v' e 

+ K +1 + <) ■ (V (2ft) - 0"(2ft))ft + - ^gp) (h ft)) (v' e+1 + v' e ) 



Finally, rearranging the sum, we obtain the final expression 

N r 



£ U ■ (V'(ft) + 40"(2ft))ft + - ft)) «i 

+ K+i + u't) ■ (^f^ - p t)) W+i + ~ £2 <+i ' (>(2ft)ft) «JVi (3-5) 
+ 2uJ ■ U"{2F 1 ){P t + ft_ x - 2ft)) ^ 

+ + «J) • (V (2ft) - 0"(2ft))ft + (^gp - ^)) (J a - ft)) (v' £+1 + v' e ) 



Similarly, the first variation of the Cauchy-Born approximation ()2.2[) is 

A ,CB, M , ^'(ft) + 2^(2ft) , , 

and the second variation is 

d 2 £ c »(y)[u,v] =ej^u' r (V'(ft) + 40"(2ft))ft + (J a - ft)) »J. (3.6) 

4. Stability of a linear chain 

In this section, we will consider a 1-D chain yp € y of atoms with nearest neighbor interatomic spacing 
Fe and study its stability with respect to two types of perturbations. First, we only consider displacements 
iieWso that yF + « € 3^; this means that the atoms can only move in the direction of the chain and yp +u 
is still a 1-D chain. These results are given in Theorems 14.11 and 14.21 and were first explicitly given in [7J . 
In the second approach, we consider 2-D displacements that allow atoms to move out of the straight line. 
These results are given in Theorems 14.31 and 14.41 We then discuss and compare the results. 

4.1. Stability of a 1-D constrained chain. Consider the 1-D atomic configuration yp S y with inter- 
atomic spacing Fe. Note that in this case ft = ft = F (see (|3.1[) for the definitions of ft and ft). We then 
have the following stability results with respect to displacements that preserve the one-dimensionality of the 
chain (cf. [7]). 

Theorem 4.1. Let yF € y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe Then 

ueu\{o} |K|||2 

Proof. This follows immediately from (|3.6[) . since for u € 14 \ {0} we have (J 2 — ft)w« = for all £, and 
therefore S 2 £ CB {y F )[u, u] = {<f>"(F) + 4 4>"{2F)) □ 

Theorem 4.2. Let yp € y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe. If<j)"{2F) < 0, then 

^g°(ttF)[".»] = r(F) + 4 0"(2F) - eW(2ft, 

uGU\{0} \\u'\\j 2 



QUASI-NONLOCAL APPROXIMATION OF LINEAR AND CIRCULAR CHAINS 



6 



where 

\\U"\\ ( 2 

ix e = inf ' ' = 2n + 0(e 2 ) as e -> 0. 

ueu\{o} \\u \\e* 

Proof. This follows from (|3.5[) . since for u e W \ {0} we have (J2 — -Pg)i4 = (^2 — -^X^+i + u 'e) = an d 
P £ + Pi-x - 2P e = for all ^, and therefore 

5 2 8*{y F )[uM = (d>"(F) + 4</>"(2F)) - e a ^"(2F)||u // ||?j. 

The result follows by applying the identity [7 . 29 

2sin(7re) 9x 

inf p = i — ^=27T + C(£ 2 ) ase^O. 

ueu\{o] \\u'\\£2 e 

□ 

Remark 4.1. VFe recall that stability means that the infima in Theorems \4-l\ and \4-2\ are positive. For a 
typical potential <f>, such as the Lennard- Jones potential, the requirement of positiveness provides an upper 
bound on the strain F for which yp is stable in the respective model. That is, it provides a bound on the 
amount of stretching the 1-D chain can undergo and remain stable. Note that the atomistic model exhibits 
slightly more stability than the Cauchy-Born approximation in the sense that the upper bound on F for which 
yp is stable is larger in the atomistic model than in the Cauchy-Born model. Asymptotically, however, as 
e — > 0, the two infima agree. Finally, we also note that the stability region is only bounded by from below, 
that is, there is technically no bound on the amount of compression the 1-D chain can undergo and become 
unstable in either model. 

4.2. Stability of a 1-D unconstrained chain. Consider again the 1-D atomic configuration yp £ y with 
nearest neighbor interatomic spacing Fe so that again F\ = F2 = F. Since under compression of the chain 
one would expect the chain to exhibit some type of buckling, we next consider displacements u £ U that allow 
the chain to become two-dimensional, and we provide stability results with respect to such displacements. 

Theorem 4.3. Let yp € y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe Then 

mf — = mm <4> (F) +4<p (2F), 

ueu\{o} 1/ ; t y j, F 

Proof. Using expression (|3.6[) for the second variation of £ CB , we immediately get 

N 



5 2 £ CB (yp){u,u}=eJ2 



e=i 



{<f>"{F) +^"(2F)) \\P e u' e \\ 2 + 0,(F)+ F 20,(2F) H(/ 2 - PKII 2 



>mmU) (F)+4<f> (2F), >\\u ||^. (4.2) 

To obtain the expression for the infimum, we will consider two types of displacements. First, for any 1-D 
displacement u £ U, we have (I2 — Pe)u'g = and 

8 2 £ CB {yp)[u,u] = ^\F) + A<t > l \2F))\\u'\\%- 

Next, consider the case of N even and the displacement u £ U such that for all £ we have u^i = Cv and 
U2i+i = —Cv for some C > and a vector v £ M 2 orthogonal to the chain yp. This case corresponds to 
creating a zig-zag deformation of the 1-D chain. In this case it is easy to see that Pgu' e = for all £ and 

6 8 (yp)[u,u} = (4.3) 

Finally, in the case of odd N, the same zig-zag deformation u defined for £ = 1, . . . , N — 1 and uat = will 
still satisfy P e u' e = for all £ and P~3|) still holds. □ 

For the stability in the atomistic model we have the following result. 
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Theorem 4.4. Let yp € y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe. If<f)'(2F) > and (f>"(2F) < 0, then, as e -> 0, 

mm\(/)"{F)+A(f>"{2F) + 0{e 2 ), tlElX if N is even, 



S 2 £ a (y F )[u,u] _ v--)-r-*v ^-j-r^ h — 

ueu\{0} |min(^"(F)+4 0"(2F) + O(£ 2 ), + 



Proof. Using expression (|3 . 5() for the second variation of £ a and observing that Pi + Pi-i — 2Pi = 0, we get 

jv 

j 2 £*(y F )[u,u]=eY, 



\<f>"(F) + 40"(2F))||P^|| 2 + ^P||(J 2 - PK|| 2 



' »(l a - P f )K +1 + ^)|| 2 - £ V(2F)||P^' + iH : 



2P 

<0"(F)+4^"(2F), 

To obtain the expression for the infimum, one can again use the displacements u and u as in the proof of 
Theorem 14.31 It is easy to verify that they satisfy (for all values of t for it, for all values of I for u if N is 
even, and for all but three values of I for u if N is odd) 

(I 2 - Pe)u e = (I 2 - Pi){u l+1 + u^-i) = 0, Piu" +1 = u'l +1 , 
P£u' £ = u e+1 +u' i = Peu" +l = 0. 

In the case of odd N, we do not have u' t+1 + u' ( = for £ = 1, N — 1, and AT, only p({^ +1 + u^) = 0, and 
therefore 

JV 

5 2 8*{y F )[u,u] = (4>"(F) + 4</>"(2P)) ||£'||| - ( £ V'(2P)) £ ]T ||P< +1 1| 2 , 



<5 2 £ a (y F )[?j,,u] = ^r-\\u'\\% if JV is even, 



<S 2 f a (y F )[u,u] = ^l||u'|| 2 2 + O(e) if jV is odd, 
F e 



and the conclusion of the theorem follows. □ 



Remark 4.2. Comparing Theorems \4-3\ and \4-4\ t° their constrained one- dimensional counterparts \4-l\ 
and \4-%\ we see that in order for yp to he stable in the two-dimensional models, additional inequalities, 



^-^ — > and — ^ - > 0, respectively, must now be satisfied in the limit as e — > 0. For typical 

potentials <f>, these inequalities provide lower bounds on the stretch F , below which the chain could undergo 
a zig-zag buckling as demonstrated in the proofs. Note that a zig-zag configuration with all nearest neighbor 
interatomic distances equal to £ and all turning angles alternatingly equal to ±27r/3 would produce a global 
minimum of the atomic energy £ a . (Following we consider the signed turning angle — n < < n 

at each atom ye defined as the angle between y' e and y'i + i, measured in the sense that negative sign of Pe 
corresponds to a clockwise turn and positive sign to a counterclockwise turn.) 

We note that in all four theorems, the inequality determining the upper bound on F is asymptotically the 
same, 4>" {F) + 4 0"(2P) > 0. Intuitively, this is reasonable to expect, since under tension, the atoms would 
tend to align themselves along a straight line, thus erasing the lowest-order difference between the atomistic 
and Cauchy-Born models, and between the one- dimensional and two-dimensional models. 

On the other hand, there is a difference between the inequalities for the lower bound on F in Theorems \4-S\ 

and \4-4\ — ^ ^ ^ > and — — - > 0, respectively. These lower bounds provide buckling thresholds 
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for the one- dimensional chains, and it is reasonable to expect them to be different in the atomistic and 
the Cauchy-Born models, since the instability mode, the zig-zag deformation, cannot be described in the 
Cauchy-Born model. Since for typical potentials <fi one has <p'(2F) > 0, the Cauchy-Born model exhibits 
more stability under compression compared to the atomistic model. The lower bound on F in the Cauchy- 
Born model corresponds to F < 1, for which also 4>'(F) < 0. If one creates a small zig-zag perturbation 
in this configuration, the nearest neighbor distance increases, thus lowering the contribution to the energy 
from the nearest neighbor interaction. The next nearest neighbor distance also increases in the Cauchy-Born 
model, resulting in an increase in the overall energy. However, in the atomistic model, the next nearest 
neighbor distance stays the same, and the overall energy is thus decreased by the zig-zag perturbation. Thus 
the requirement of<f)'(F) > in the atomistic model seems quite reasonable. 

5. Stability of a circular chain 

In this section, we will consider uniform circular configurations yF € y oi N atoms with nearest neighbor 
interatomic spacing Fe and study their stability with respect to 2-D perturbations u £ U. The atoms lie on 
a circle of radius R that satisfies 

Fe = 2i?sin(7re). 

The distance between next nearest neighbors of the chain y is 2Fecos(7re). Therefore, we have 

Fi = F and F 2 = F cos — — F cos (ire). 

N 

We note that due to the symmetry of the circle the only forces on the atoms in both the atomistic and the 
Cauchy-Born model are in the radial direction, i.e., the direction normal to the circle. It is straightforward 
to obtain from equations ()2.1|) and (|2.2|) that these forces vanish in the Cauchy-Born model if 

<j>'(F)+2<j>'(2F)=0, 

and they vanish in the atomistic model if 

<^(Fi) + 2cos(7r£)^'(2F 2 ) = or <j)'{F) + 2 cos(Tre) </>'(2 cos(tts)F) = 0. 

We remark that the above equations arc special cases of the vanishing of the first variations of the respective 
energies, S£(y F ) = 0, as given in Section |3l For Lennard- Jones type potentials <fi, these equations always 
have a solution, and thus both the atomistic and the Cauchy-Born model possess a circular equilibrium, 
although in general with different radii R a and R CB . Using the inverse function theorem, the difference 
between the radii can be seen to be R a — R CB = 0(e 2 ) as e — > 0. 

We now have the following result for the stability in the Cauchy-Born model. 

Theorem 5.1. Let yp G y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. Then, as e — > 0, 



inf 



S 2 £ CB (y F )[u, 



u. 



min i 4>"{F) + 4 <t>"(2F), + J K if N is even, 



" ew \<°> " lmin(^(F)+4^(2F), ^ +^'^ +0(e) | if N is odd. 

Proof. Using expression (|3.6j) for the second variation of £ GB , we immediately get 



N 



S 2 £ CB (y F )[u,u]=eJ2 



( = i 



(<ro + 4 0"(2F)) np^ii 2 + fin+i&i m _ Pl)u > M * 



4>'{F) + 24>'(2F)\. 



>min|0"(F)+4 0"(2F), yv ' ' p Y ' J *f\\u'\\%- (5.1) 

To obtain the expression for the infimum, we will again consider two types of displacements. First, if 
U{ = Cye for some C > and for all £, a case that corresponds to the pure expansion of the circle, then 
u eW, (I 2 — Pe)u' p = 0, and 

S 2 £ CB (y F )[u,u]^(cf>"(F)+4^(2F))\\u'\\l. 
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Next, consider the case of N even and the displacement u such that ii 2 £ = Cy 2 £ and u 2 t+i = —Cyit+i for 
some C > and for all I. This case corresponds to creating a zig-zag deformation of the circle. In this case 
it is easy to see that u eW, Pgu' e = 0, and 

d£ (y F )[u,uJ = IMI^- (5.2) 

r £ 

Finally, in the case of odd N, the same zig-zag deformation u defined for £ = 1, . . . , N — 1 and ujv defined so 
that X^=i = will have two segments, one joining the first and JV-th atoms and one joining the (N — l)-st 
and iV-th atoms, for which Peu' e ^ 0, thus creating a 0(e) perturbation in (|5.2[) . □ 

Before we address the stability in the atomistic model, we note that for circular arrangements P( + Pi_i — 
2Pt ^ for any I. However, it is easy to check that due to the geometry of the circle, we have 

|| (P e + Pz-i - 2Pt)w\\ = 2\\w\\ sin 2 = 2||io|| sin 2 (ire) < 2tt 2 £ 2 |HI for all w £ R 2 . (5.3) 

We now have the following stability result. 

Theorem 5.2. Let yp £ y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. If <j) £ C 3 (0, oo), (f>'{2F) > 0, and 4>"(2F) < 0, then, as e — > 0, 



inf S 2 £ a {y F )[u,u 



min i cf>"{F) + 4 0"(2F), < ^p- \ + 0{e 2 ) if N is even, 



uew\{0} |min|0"(F) +40"(2F), ^p- +0(e)J + C(e 2 ) if N is odd. 

e 2 

Proof. First, recall the definitions F\ = F and F2 = Fcos(7re), so \Fi — F\ < — . Using the smoothness of 
4>, there exists a constant C§ > 0, independent of e, such that 



b"(2F 2 )-<f>"(2F)\, 



c/>'{2F 2 ) (j>'{2F) 



2F 2 2F 

Using expression (|3.5[) for the second variation of £ a , we now have 



<C^e l . (5.4) 



N 



5 2 £ a (y F )[u,u]=eJ2 



(4>"(F) + 40"(2F))||P^|| 2 + ^p-\\(h P e )u' e \\ 2 

F) -\\(I 2 P e )(u> +1 + u>)\\ 2 - e 2 ^(2F)\\P e u>; +1 \\ 2 



2F 

+ 2u' r U"{2F)(Pt + Pe-i- 2Pt) 
+ {<f>"{2F 2 )-<t>"{2F)\\P l {u' i+1 +u' l )\ 



'M|2 



b'{2F 2 ) cj>'{2F) \ 
2F 2 2F J 



\\(I 2 -PlM +1 +u' e )f 



Using the assumptions <p'(2F) > and <f>"(2F) < 0, applying (|5.3[) . and (|5.4[) together with the triangle 
inequality, we get 

,m Win.. „,2 



<5 2 £ a (y F )[ U , U ] > min + 4 0"(2F), ^ j 

- 4^ 2 e 2 |0"(2F)| -4 £ 2 <7 ||</|||. (5.5) 

To obtain the expression for the infimum, one can again use the displacements u and u as in the proof of 
Theorem 15. II It is easy to verify that they satisfy 

P t u'l +1 = 0, (I 2 - Pt){u' t+1 + ^_i) = 0, 

P t u'[ +l = 0, (I 2 - Pi)(u' l+1 + u 1 ^) = 0, 

except for two values of ug when N is odd, and the conclusion of the theorem follows. □ 
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Remark 5.1. Asymptotically, as e — > 0, the results in Theorems I5.il and \5.S\ are equivalent to those in 
Theorems \4-3\ and \4-4\ The requirement of stability gives two inequalities, and these inequalities produce an 
upper and lower bound of F for typical potentials <j). Yet again, we can interpret the upper bound as a limit 
on the amount of stretching the chain can undergo before fracturing, while the lower bound can be interpreted 
as a limit on the amount of compression before the chain starts buckling and creating zig-zag segments. 



6. Stabilization by a bond-angle energy 

Many-body empirical potentials generally include the effect of bond angle in addition to two-body inter- 
actions [3J|3n]. We will study a simple bond-angle energy in this section and show that it suppresses the 
buckling modes under compression when added to the atomistic energy (|2.1[) and the Cauchy-Born energy 



(6.1) 



Wc will consider a bond-angle energy of the form 

£ h (y) =e$>(i -cos&), 



where a > is a constant and — 7r < f3i < it is the signed turning angle of the chain y at the atom ye. Recall 
from Remark 14.21 that /3t is the angle between y' g and y'g +1 and that negative sign of Pi corresponds to a 
clockwise turn and positive sign to a counterclockwise turn. We note that 



cos (3i 



Ve+i 
\y'e+x\ 



Well 



To see the effect of this term on the overall energy, we first compute the first and second variations of £ b . 
For the first variation, we have the following result. 

N r / /l l 



5£ h {y)[u] = aeY, 

N r 

E 

N 

E 



\vi+i\ 



Vi+i 



Vt 





( u 'i y 




y'i+i 




\\WiW h 


' t \\J HWJ 


\\y'e + i\\\ 



ae 



cos (3 t 



( u £+i y'e+i 



(6.2) 



ae 



UwJ+ill W t+ i 

(cos/3 f + cos/^_i) 



\Vi\\ WtW 



\Ve+i\ 



vtW WtW 



y'e-i 

WvLi 



\y'e+i\ 



where in the last step we rearranged the sum using the periodicity of the chain. We now notice that if all 
turning angles are the same, that is, if fii = /3 for all t, then we also have 



and 



5£ h (y)[u]=aeJ2 



Vt+i 
\Vt+l\ 



2cos/3 



Vl-i 



= 2cos/3/^, 

WvtW 



(6.3) 



WA ' ll^l 



2cos/3 



0. 



This means that, among others, straight chains and uniform circular chains are critical points of this bond- 
angle energy. 

To compute the second variation of £ h , we start with the first expression from (|6.2j) and use the product 
rule twice. After some simplifications we obtain 
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Ve+i _ Vi Vi _ y^+i 



ll^ll ll^ll lltti+ili; ll^ll 



IM+J V ll^+ill ll^+ili; llai+ill 



mii bin; w 

The first two and the last two terms can be combined due to the periodicity of the chain, and we obtain 



6 2 £ h (y)[u,v] =ae^ 



' y' e+1 y'i.x \ y' t y' e ( y' e+1 yl 



HW WWv'i+A WvUW) \\v\W m\ VIM+ill WvUW) ) W t \ 

+ (cos/3, + coflft_i) i^fji • (h - 3 ,4n ® FTii) 

Ikll V Ikll Ik- 11/ 



Mil v iiv2ii w m 

Finally, for a deformation yp for which all nearest neighbor interatomic distances are Fe and all turning 
angles are /3, we can again use (|6.3[) to get 



^ b (y F )[ M ,,] = e fv« • fe?£(I a - ft)) v>. 



/ y rv rns /-I \ 

(6.4) 



Note that if the turning angles satisfy — tt/2 < (3 < x/2, then 

S 2 S h (y F )[u,u] = ^P^||(I 2 - Pi)u' t \\% > 0, (6.5) 

and if < tt/2, then 

(5 2 £ b (2/ F )[u, u] = if and only if (7 2 - ftK = for all 

which corresponds to pure expansion or compression of the circular chain yp. Clearly, in this case the angles 
j3 do not change and the bond-angle energy does not either. 

Using the second variation (|6.4j) of the bond- angle energy £ h , we can now re-state Theorems l4.3[ |4~41 15.11 
and 15.21 for the augmented energies £ CB (y) + £ h (y) and £ a (y) + £ h (y). 

Theorem 6.1. Let £ a ' b {y) = £ a (y) + £ b (y) and £ CB ' b (y) = £ CB (y) + £ b (y). Let y F G 3> denote the 1-D 
configuration of atoms with nearest neighbor interatomic spacing Fe so that all turning angles are 0. Then 

inf ^^^(ffi|lf( 2 F) nF) + 2<f>> m , 2a 



ueu\{oy \\u'\\% r ^ F 2 V 

e 

and if </)'(2F) > and <p"(2F) < 0, then, as e -> 0, 

5 2 5-- 6 ( I , J ,)[«,tt] Uin|0"(F)+40"(2F)+O( £ 2 ), ^ + ^| z/TVzse^, 
«e«\{0} | min + 4 ^ (2F) + 0(e 2 ); ffi) + |a + ^ j , /jV fa ^ 

Lei y F & y denote the uniform circular configuration of N atoms with nearest neighbor interatomic spacing 
Fe so that all turning angles are f3 e = 2it/N = 2ire. Then, as e — > 0, 



inf 



1 min (ffl +4f-(2F), ±1/03. + 2 °^ A + Q( e ) } if N is odd, 
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and if(f> € C 3 (0,oo), <f>'(2F) > 0, and <p"(2F) < 0, then 



inf 



S 2 £^(y F )[u,u] Jmin|^(F)+4^(2F), ^ + ^p^}+0(e 2 ) i/7V e.en, 

RE = 1 ,^ in |^ (jF) + 4 ^ (2F); + + G(e) J + 0(£2) 4/ jv & oM 

Proof. The claims follow immediately from Theorems 14. 31 [4~4l 15.11 and 15 .21 and the second variation (|6.4p of 
the bond- angle energy £ b . □ 

. . +2 0'(2F) 

Remark 6.1. Note that in Theorem \6 . 1\ the expressions — in the formulas for the infima of 

the Cauchy-Born energies and the expressions — ^ - in the formulas for the infima of the atomistic energies 
2(X COS fj£ 

are augmented by — > if all turning angles are the same and satisfy \/3 e \ < n/2. Also note that, for 

typical potentials such that (j)'(F) <0if0<F<l, the additional bond-angle term lowers the lower bound 
on F below which the chains can undergo the zig-zag buckling. Note, however, that for the Lennard-Jones 

, , <p'(F) ' , , , 2a cos ^ e 
or Morse potentials, the term — — — is negative and dominates the term as F — > 0, so no matter 

how large a is, there is always a positive lower bound for the region of stability of the deformation y F . 

Remark 6.2. The bond-angle energy (|6.1[) given by £ b (y) = ey]^ a(l — cos/3^) is meant to model in 
one- dimensional chains the resistance to transverse displacement in two-dimensional graphene and carbon 
nanotubes. There is an energy cost to transverse displacement in tri-bonded graphene and carbon nanotubes 
when modeled by popular potentials such as \30\j which penalize the bond-angle deviation from 2it/3. 



7. Modeling errors 

In this section, we will consider applying an external periodic load / € U to a deformation y F € y and 
study the error of the Cauchy-Born model. More specifically, let u a € U and u CB S U solve the linearized 
equations 

£ a (VF)[v] + 5 2 £ a {y F )[u a ,v] = (f,v) for all v e U, 
S£ CB {y F )[v}+5 2 £ CB (y F )[u CB ,v} = (f,v) for all v G U. 
We define the modeling error of the Cauchy-Born approximation, r, via the duality relationship 

(t,v) : = 6£ GB (y F ) [v] + 6 2 £ CB (y F ) [u 8 , t>] - </, v) 

= {S£ CB (y F )[v] - S£ a (y F )[v]} + {5 2 £ CB {y F )[u\v\- 8 2 £ a (y F )[u a ,v]} for all veU. 

Since the solution u will play no role in the analysis of the modeling error, to simplify the notation in the 
rest of this section, we will sometimes suppress the superscript and simply write u instead of u a . However, 
in the statements of the theorems, the proper notation will be used. 
We compute that we have for all v € U 



"(2*1) 



2Fi 2F 2 I v a ^ l+L > 2F X 



S£ CB (y F )[v]~6£ a (y F )[v]=eY, 

and thus observe that 

5£ (y F ) — 5£ & {y F ) = if y F is a linear chain since Fx = F 2 and y'p = 0, 

but that 

S£ CB (y F ) — S£ a (y F ) ^ if y F is a circular chain. 



v' e (7.2) 
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7.1. Modeling errors for linear chains. Using the second variations of £ a and £ CB given in (|3.3[) and 
(|3.6p . respectively, we now have since S£ (yp) — 8£ & (yp) = for linear chains that 



(t,v) = 5 2 £ L,u (y F )[u, v] - 5 2 £%y F )[u,v 

N r 

= e E 



Pi 



^•^"(F 1 )P, + ^l(/ 2 -P,))^ 



If we define the matrices 



A t := 4>"(2F 1 )P, 



(uj+i+uj)- U"(2F 2 )P< + 



6'(2Fi) 



^2P 2 ), 
2P 2 



(J 2 -P<) W+i+^) 



(7 2 - P>) and if 



(2P 2)j P, + ^)(J 2 -P,; 



2i?i V w ...... 2^ 

cancel the terms on the second line and use identity (|3.2[) . after rearranging the sum we get 



N 

(t,v) =s^2 
t=l 

N r 



4 u' e ■ A e v' e - 2 u' e+1 A e v' e+1 - 2 u' e A e v' e + e 2 u" +1 Aiv'[ +1 



E 



e 2 u" +1 ■ A t v'l +l -2v! r [A l +A l _ 1 -2AA v' t 



We can now sum by parts and simplify to get 

N 



{T, V) = -£ 



E 



e 2 A 



eu e+1 



A e - A t _ x 



2 [At + A t -i - 2AA u' e 



(7.3) 



(7.4) 



Let us now consider a 1-D atomic configuration yp £ y with interatomic spacing Fe. In this case we 
have F\ = P 2 = F, and also all of the projection operators Pi and Pi are the same for all £, so we can write 
Pi = Pi = P and Ai ~ Ai = A. The expression (|7.4p for the modeling error then simplifies to 



and we immediately have the following theorems, one for the 1-D constrained chain and one for the 1-D 
unconstrained chain. 

Theorem 7.1. Let yp £ y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe, let f £ U, and let u a <EU satisfy 5£ a {yp)[v\ + S 2 £ a (yp)[u a , v] = (/, v) for all v £ U. The modeling error 
of the Cauchy-Born approximation, t, then satisfies the inequality 

\\r\U<e 2 \ ( t ) "{2F)\\\{uT\\q- (7.5) 

Theorem 7.2. Let yp £ y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe, let f <EU, and let u a £ U satisfy S£ a (yp)[v] + S 2 £ a (yp)[u a , v] = (/, v) for all v £ U. The modeling error 
of the Cauchy-Born approximation, t, then satisfies the inequality 



|t||. < £ 2 max||0"(2P)|, 



4>>(2F) 



2F 



l(«TII< 



(7.6) 
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7.2. Modeling error for circular chains. The situation is different in the case of a uniform circular chain 
Uf £ y with interatomic spacing Fe. In this case, S£ (uf) — SE^d/p) ^ 0, but recalling the constant Ca, 
from (|5.4p . we obtain from (|7.2I) 



5£^{y F )[v]-5£*{y F )[v] 



N 

■E 



b'(2F 1 ) ( j>'(2F 2 ] 



(4^ 



< 



2Fi 2F 2 
e 2 C,(4||^|| £ ,+e 2 ||^'||^) 

4C (l + ^ 2 e 2 ) J F + 2^ 2 |0'(2F)| 



yF,£+i 



«£'(2F) 




2F 


Il2/£'lk 



2Fj 



= e 2 c K !| v '|| e 



(7.7) 



where 



(7.1 



C k =4C^I + tt 2 s 2 )F + 2tt 2 |0'(2F)| 
since = F and = 4Fe~ 2 sm 2 (7re) < 4F?r 2 . 

For the modeling error term 5 2 £ CB (uf)[u, v] — 5 2 £ a (yF)[u, v], there is no cancellation in (|7.4|) as in the case 
of the straight chain, but the expressions Ai — Ai_\ and Ag + Ai_\ — 2Ai are of order e, and e 2 , respectively. 
Therefore, the modeling error of the Cauchy-Born approximation is again of order e 2 . 

Theorem 7.3. Let yp £ y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe, let f <EU, and let u a eW satisfy 5£ a (yF)[v] + S 2 £ a (yF)[u a , v] = (/, v) for all uetf. 
//</>£ C 3 (0, co), then the modeling error of the Cauchy-Born approximation, t, satisfies the inequality 

\\r\U < C K e 2 + (C l£ 2 + C 2 e A ) (\\(u a )"'hl + + IKO'll*?) 

where 



d = max 1 1 <j>"{2F) |, 



<f/(2F) 



2F 



, 4tH 



\2F) 



b'{2F) 



2F 



12 C d 



C 2 =4ttC , 



(7.9) 
(7.10) 



is the Lipschitz constant defined in (|5.4I) . and C K is defined in 

Proof. Using the definition (|T. 3[) of Ag and Ai, we have, after some manipulations, 

V{2F 2 ) 



At 



2Fo 



Pi — Pe-i , 



and 



Ae 



Ai-i - 2A t 



'(2F, 



</>'(2F 1 ) 



2F 1 

"(2F 2 )-<f>"(2F 1 )) 



4>'{2F 2 ) (f>'{2F l y 



2F 2 



2Fi 



Pi + Pt-i 



It can be easily verified that for any vector w <E R we have (cf. ([5 

2ir 

\\(P t - P e -i)w\\ = \\w\\sm— = ||w||sin27T£ < 27rc||w)||, 
\\(P e + Pe-! - 2P e )w\\ = 2\\w\\ sin 2 j- = 2\\w\\ sin 2 (tte) < 2ir 2 e 2 \\w\\. 

Recall that Fi = F and F 2 = F cos (7re), and from (|5.4p that there exists a constant > independent of 
e such that 

4>'{2F 2 ) cj ) '(2F 1 ) 



max||0"(2F 2 )-0"(2F 1 )|, 
Therefore, from (|7.4[) and (|7.7|) we obtain 



2F, 



2Fi 



Ml* <e 2 max \\(j>" {2F 2 )\, 



V{2F 2 ) 
2F 2 



\\{uy\U+2ne 2 



b"(2F 2 



ft{2F 2 ) 
2F 2 



\(u a )"\V 
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2ir 2 e 2 



"(2F 1 ) 



b'(2F 1 ) 



2F t 
f(2F 1 ) 



+ 6C^ ||(u a )'ll« + C K e 2 



2n 



4tt z 



a" 



(2*i) 



'(2Fi 



2F 1 



2Fi 
0'(2F X ) 



2Fi 



- £ 2 

4tt Q e 2 
f 12 Q 



(«Tlk 



and the statement of the theorem follows. 



□ 



Remark 7.1. The above second-order modeling error estimates 0{e 2 ) for the Cauchy-Born approximation 
were the result of the symmetric treatment of the interactions between next nearest neighbors. A more explicit 
treatment of the second-order modeling error for the Cauchy-Born approximation for linear chains was given 
in JEjj. 

We note, however, that the second-order modeling error estimates C(e 2 ) for the Cauchy-Born approxima- 
tion require that \\(u a )'"\\(2 be bounded uniformly ine, which is not the case for the approximation of atomistic 
configurations with defects. This lack of accuracy of the Cauchy-Born approximation for problems with de- 
fects is the motivation for the development of atomistic-to- continuum methods such as the quasi-nonlocal 
method which attain 0(e 3 ^ 2 ) accuracy for problems with defects JB,[1~4 



8. Error analysis for the Cauchy-Born approximation 
From the definition (|7.1[) of the modeling error of the Cauchy-Born approximation, we have 



(r,v) = 8£ (y F )[v] + 8 £ (y F )[u a ,v] - (f,v) 



5£ (Vf)[v] + 8 £ (y F )[u a ,v] — 8£ [y F )[v] — 8 £ (y F )[u 



,CB 



= 8 2 £^{y F )[u\v]-5 2 £^(y F )[ 



,CB 



u CB ,v] 



Setting v 



,CB 



= 5 2 £^{y F )[r 
in (|8.1[) above, we obtain that 



for all neW. 



8*£^{y F )[u a 



u CB ,u* 



,CBi 



(r, u a - u j 



(8.1) 



■2) 



Let us also define the Cauchy-Born stability constants (see Theorems 14.11 14.31 and 15. ip 

7i = <f>"(F) +4(f)"(2F) and 



l2 = mta {,» OT+ 4*»(2 n «a±^}. 



Using the stability results and the expressions for the modeling errors, we now have the following theorems 
for the errors of the Cauchy-Born approximations. 

Theorem 8.1. Let y F € y denote the TD configuration of atoms with nearest neighbor interatomic spac- 
ing Fe. Given f G U, let u a ,u CB € U satisfy 5£ a (y F )[v] + 8 2 £ a (y F )[u a ,v] = (f,v) and 8£ CB (y F )[v] + 
5 2 £ CB (y F )[u CB ,v] = (f,v), respectively, for all v £ U. If "fi > 0, then 



,CB\) 



7i 



(u a )'"\\ 



Proof. Combining (|8.2j) with (|4.1|) and fjT. 5[) . the theorem follows from the inequality 



7i||K - O ||*2 < 5 z £^{y F )[u* - u 
' = (r,u a -u CB ) 

< ||t|U||(«--« ct ) / || 



,cb 



,CBi 



< e 2 i0"(2F)i!iKrii£ ? ilK-« CB )'ll^ 
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□ 

Theorem 8.2. Let y F £ y denote the 1-D configuration of atoms with nearest neighbor interatomic spac- 
ing Fe. Given f g U, let u a ,u CB g U satisfy S£ a (y F )[v] + 8 2 £ a {y F )[u a ,v] = (/,«) and 6£ CB (y F )[v] + 
5 2 £ CB (y F )[u CB ,v] = (f,v), respectively, for all v g U. If 72 > 0, then 

b'{2F) 



e 2 max<j \<j>"{2F)\, 
\\{u a -u CB )'\\ P < 



2F 



\(u a Y"\ 



72 

Proof. Combining (|8.2p with (|4.2[) and (|7.6p . the theorem follows from the inequality 

CBwi|2 ^ r2cCB/„, u„a „,CB „,a „,CBi 



72 ||(u a - ^ B )'|| 2 2 < 5 2 £^ B {y F )[u* - u™,u 



(T,U a ~U CB ) 

< ||r|U||( U a -,i CB 
max\\(j)"(2F)\, 



cj>'{2F) 



2F 



||(«T'llell(« B -«^yil/ 



□ 



Theorem 8.3. Let yp g y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. Given f g U, let u a ,u CB g satisfy 5£ a (y F )[v] + S 2 £ a (y F )[u a , v] = (f,v) and 
5£ CB {y F )[v] + S 2 £ CB (y F )[u CB , v] = (/,«), respectively, for allvEU. If <f> £ C 3 (0, 00) and 72 > 0, ttera 

IKu'-u 03 )'^ < 72" 1 [C K e 2 + (C l£ 2 +C 2£ 4 ) (||(«Tlk + \\(uT\\ ei + \\(u a Y\\ ei )} , 

where C\ and C2 are defined in (|7.10l) and C K in (|7.8I) . 

Proof. Combining (|BT2j) with (f5TT]l and ([715]) . we obtain 

72 ||( W a - W CB )'||| < 6 2 £ CB (y F )[u* - U CB ,u a - U CB ] 
= ( T ,u a -u CB ) 
<||r|U||(«»-« GB )^ 

< [ev 2 + (cV + c 2 e 4 ) + + ll(« a )'lk)] ll(^ a - « CB )'lk, 

and the theorem follows. □ 

Remark 8.1. We note that the error \\{u a — u CB )' \\p will generally not be small for compressive F < 1 such 
that the atomistic model is not stable ((f) 1 '(F) < 0), but the Cauchy-Born model is stable ((£>' (F)+2(j)' (2F) > 0), 
because || (u a )"'||^2 can be expected to be large. 

9. Error analysis of the models including the bond-angle energy 



The error analysis for the case when the bond-angle energy £ b as defined in (|6.ip is included in the models 
is straightforward. The stability results arc given in Thcorcm l6.ll Since we only need an estimate from below 
on the Cauchy-Born stability constant, we can use again 72 as defined in Section [51 

There is no contribution of the bond-angle energy £ h to the modeling error, t, since t is given in (|7.ip by 

(t,v) = {S£ CB ' h (y F )[v} — S£ a ' h (y F )[v}} + {S 2 £ CB ' b (y F )K> b , v] - 5 2 £^{y F )[u^, v}} 

= {6£ CB (y F )[v\ - 5£*(y F )[v]} + {5 2 £ CB (y F )[u a > h ,v] - <5 2 £ a (y F )K b , v]} , 

where 6£*' h (y F )[v] + S 2 £ li - h {y F )[u^ h ,v} = (f,v) and 5£ CB ' h {y F )[v] + S 2 £ CB - h (y F )[u CB ' h , v] = (f,v) for all 
v g U. Therefore, the estimates (|7.5[) , (|7.6p . and (|7.9p for the modeling error still hold with u a replaced by 
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Wc therefore have the following theorems, which are analogous to Theorems 18.21 and 



Theorem 9.1. Let y F G y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe. Given feU, let u a >\u CB - b e U satisfy S£ a - b {y F )[v] + S 2 S a - b {y F )[u a -\v] = (f,v) and 5£ CB > b (y F )[v] + 
S 2 £ CB ' b (y F )[u CB ' h ,v] = (f,v), respectively, for all v E U. If 72 > 0, then 



II (« 



a. b 



. CB,b\t 



y\k < 



£ 2 max W'(2F)\, 



b'{2F) 



2F 



72 



w(u a - b r\\ ei . 



Theorem 9.2. Let y F € y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. Given f E U , let u a ' b ,u CB ' b E U satisfy 5£°" b (y F )[v] + 6 2 £ a ' b (y F )[u a ' b , v] = (f,v) 

V 

then 



and 6£ CB - b (y F )[v] + S 2 £ CB > b (y F )[u CB - b ,v] = (f,v), respectively, for allvEU. If <p E C 3 (0, 00) and l2 > 0, 



]{u a,b „ u CB, b y hi < 7 -l [Ck£ 2 + {Ci£ 2 + ^^4) (||( u «,6)'"|| + \\( U «>»)"\\ ei + || (tl^)'^)] 



where C\ and C2 are defined in (|T. 10[) and C K in (|7.8 



10. THE QUASI-NONLOCAL APPROXIMATION 

In this section, we consider a combination of the atomistic approach in one region and the Cauchy-Born 
approximation in its complement. Specifically, for an integer 1 < K < N, we define the set 21 of nonlocal 
atoms ( "nonlocal" or "atomistic" region) and the set £ of local atoms ( "local" or "continuum" region) as 

21 = {1, ... , K} and € = {K + 1, . . . , N}, 

extended by periodicity. 

For y E y, we define the quasi- nonlocal energy, £ ^ NL (y), as a combination of the atomistic energy in the 
nonlocal region and the Cauchy-Born approximation in the local region: 



£ QNL ( y ) = £ £ 



«G2l 
N 



Wi\\) + We+i + y' e \ 



E 



1 



<KH\\) + o (0(2M) + 0(2||^ +1 ||)) 



N 



= eJ2^\\y' e \\) + sJ2H\\y' e+ i + y'A) + ^ E U^U\\) + 4>(ny'i + x\\))- 

1=1 i=\ l=K+l 

We can rearrange the sums to resemble the atomistic and Cauchy-Born energies in the following way: 

N K N 



£^ L ( y ) = eJ2 <KHW) + £ E MM+i + «<H) + e E # 2 IMID + \ ^ 2 IMH) + \ ^IbWilD- 



1=1 



l=K+2 



When studying the quasi-nonlocal method and its properties, the following seminorms of displacements 
u EU will be used: 



K 



1/2 



IKIkw= e EKI 



. \\ u hue) 



(e £ H\A . \\u'\\ em =(e\\u[\\ 2 +e\\u> K+1 \\A 
\ e=K+2 / ^ ' 



1/2 



Note that since the primes denote backward differences, 1 1 it' 1 1 ^.2 (21) i s a seminorm over the bonds between 
nonlocal atoms, ||u'||^|(c) is a seminorm over the bonds between local atoms, and ||u'||^2p) is a seminorm 
over the interfacial bonds. Also note that 



l«'ll? a = Ik" 2 



«(2l) 



l«'l|? 2(e) + ll«"' : 



QUASI-NONLOCAL APPROXIMATION OF LINEAR AND CIRCULAR CHAINS 18 

10.1. Variations. The first variation, 5£ < ~ tNh (y)[u], is 

SS^(y)[u) =ejt^m y > . u > + + y>) • K +1 + u' e ) 

and the second variation is 

+ ef>J +1 +u' e ) ■ U(H+i+vl\\)Pt + M+i +v't) 

p — I V WVi+i "+■ % II / 



«=i 

JV 



=A+2 



£ ^■(40"(2||^||)P, + 2^M(/-^))^ 



+ (4 *"(2||yi ||)Pi + 2 - A)) «i 

+ f <+i • (4^(21114+111)^+1 + 2 ^ii^j 1 '' V ~ ftf+i)) v' K+l . 

Assuming again that the strains F\ and F2 defined in (|3.1[i are independent of £, we can rewrite the second 
variation as 

S 2 £^(y)[u,v] = e5>i • (V(ft)ft + ^(I-P*)) «J 



=1 

K 



e^K+i + «i) ■ ( <^"(2ft)ft + ^% - ft) ) W+i + «i) 



e X ^•(4f(2F 1 )Pi + 2^(/-P < ))^ 

£=A+2 



+ £ «i • (4 ^(2Fi)P! + 2 - Pi)) v[ 

+ \ <4+i ' U<f>"(2F 1 )P K+1 + 2 - Pa+i)) v' k+1 . 

If we now apply the same manipulations preceding the derivation of (|3.5[) . rearrange the sums suitably, and 
collect interface terms, we obtain 

S*£^(y)[u,v] = £ f>£. (V'(Pi) +4^(2^))^ + - P t )) v[ 

+ eYl[« + i + «i) • (^P^ ~ ^)) W+i + 4) " ' (V(2ft)ft) <+i 

A 

+ 2eX"i ' (0"(2Pi)(ft +ft-i - 2P £ )) i£ 



QUASI-NONLOCAL APPROXIMATION OF LINEAR AND CIRCULAR CHAINS 



19 



K 

1=1 



J2 W+i+«;)-((*"( 2 - F i)-<*"( 2f i))«) (v't+i+v't) (10.1) 

1 + K tl+ «;,.((^-*|^)(/ 2 -«)) Wtl+ »;, 



e 

e=K+2 



2^ Uf[[<P i F i)+^9 {2F 1 ))P i + — (h-Pi 



+ eu[. (yW + ^"{2F 1 ))P 1 + M^ffl (J 2 - ^ 

+ ^Wi ' ((0"(F 1 ) + 4 0"(2F 1 ))Pi f+ i + ffiO+ffi^ _ Pjf+l) ) 

+ 2e U ; • (^"(2F 1 )(P 1 - A)) «i + 2eu' K+1 ■ (> \2F 1 ){P K+1 - P K+1 )) v' K+l . 

Thus we see that the second variation can be roughly decomposed into the second variation of the atomistic 
energy (first five lines in (| 10 . 1 p ; cf. (|3.5p ) and of the Cauchy-Born energy (the sixth line; cf. (|3.6[1 ), with 
some terms arising at the interfaces between the atomistic and continuum regions (last three lines). 

10.2. Stability of the quasi-nonlocal approximation. Using the second variation of <5 2 £Q NL given in 
(|10.1[) , we can now give sharp stability results similar to those in Sections 2] and [5] for the atomistic and 
Cauchy-Born models. 

For the 1-D chain, yp £ y, we have that p> = p> = P for all I and F\ = Fi = F. In addition, for the 
constrained chain, for which the displacements satisfy u €lA : we also have (I — P)ui — for all i. Therefore, 
we immediately have the following stability result for the 1-D constrained chain. This result is consistent 
with that derived in [7]. 

Theorem 10.1. Let yp £ y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe. If<j)"(2F) < 0, then 

ueU\{0} II" lua 

Proof. Using expression (|10.1j) for the second variation of £^ NL , we have 

N K 

^ N ^^)[ U , W ]= £ ^(^(P)+40''(2P))||<|| 2 - £ ^ e V'(2P)||< +1 || 2 



=i 



>(<//'(P)+40"(2P))||u'|||. 

The sharp estimate for the infimum follows by choosing a test function u € IA supported away from the 
atomistic region. □ 

For the 1-D unconstrained chain, we still have Pi = Pi = P for all £ and Pi = F-2 — P, but since the 
displacements are no longer one-dimensional, we no longer have (/ — P)ui = for all i. The stability result 
for the 1-D unconstrained chain follows. 

Theorem 10.2. Let yp € y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe. If<f)'(2F) > and <p"{2F) < 0, then 

8 2 £Q NL {y F )[uM > min|^(P)+4^(2P),ffl|||u'||| (gt) 

+ mm.l<j> (P) + 40 (2P), >IMI^(e) 



+ mi .{r(F) + 4^(2 f ) /'' F ) + /P F ) }||„-|||, 3 , 



QUASI-NONLOCAL APPROXIMATION OF LINEAR AND CIRCULAR CHAINS 



20 



> min^ 4>"{F) + 4 0"(2F), 



<//(F) 



for all u G hi . In addition, as e — > 0, 

5 2 £Q NL {y F )[u 1 u] 



inf 

«GW\{0} 



\\u'\\ 



l\\2_ 



^"(F)+4 0"(2F),^ + O(e) 



Proof. Using expression ([10. lj) for the second variation of £Q NL , we have 



K 



1=2 



\F)+Acp"{2F))\\Pu' l \\ 



4>'{F) 
F 



-PH\\ 



K 



2F 



\\(I 2 -P)(u' e+1 +u' e )\ 



/\||2 Jliii 



<j>"{2F)\\Pu'l +1 \ 



l=K+2 



(0"(F) + 40"(2F))||P<|| 2 + ^ (F)+ F 2 ^ (2F) ||(J 2 - PKf 



■ £ 



+ £ 



{<t>"{F) + A<i>"{2F))\\Pu' i r + 
{ ( j ) "{F)+A^{2F))\\Pu' K+1 \\- 



</>'(F) + 2</>'(2F) 
F 1 



q/(F) + 2^(2F) 
F 1 



(h-P)u' K+1 \ 



Due to the assumptions 4>'(2F) > and (f>"(2F) < 0, we can drop the second sum and estimate the remaining 
terms to get 



S 2 £^(y F )[u 7 u] > wiaU'iF) + 40"(2F), ^U||u'|| 



F 



\%{%) 



7' 



+ min<U"(F) +4 0"(2F), 



c//(F)+c//(2F) 
F 



,'l|2 



1^(3) 



Since </>'(2F) > 0, the smallest of the minima above is the first one, and the first part of the theorem follows. 

To obtain the sharp result for the infimum, consider again the displacements ii G U and u G IA as in the 
proofs of Theorems 14.31 and 14.41 If the support of u G U is chosen so that = for I = 1, . . . , K, then we 
obtain 

S 2 £^ h (y F )[u,u] = (c/>"(F) +40"(2F))||u'|||. 

On the other hand, if the zig-zag displacement u G U is chosen so that its support is contained in the 
atomistic region (atoms 2, ...,K — 1 to be specific), then there is no contribution to <5 2 £^ NL from the 
continuum or interfacial regions, however, there is a 0(e) contribution from the and u' e+1 + u' e terms 
near the boundary of the support. We then have 



and the theorem follows. 



□ 



Finally, in the case of the uniform circular chain, none of the projection operators are the same, nor 
Ft = Fz- However, using the techniques used in Section [5j we have the following result. 

Theorem 10.3. Let y F G y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. If 4>'(2F) > and <j)"(2F) < 0, then 



5 2 £^ NL (y F )[u, u] > min L" (F) + 4 0"(2F), ^1 



,'l|2 



Uf(3t) 
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+ minU"(F) +4 0"(2F), 



+ mm^"(F)+4:cf,"(2F), 

-£ 2 (47r 2 +4C )|K||! (2l) 
-s(27r\<j>"(2F)\ + 2eC 4 ,)\\u'\\% c( i ) 

> m m|^(F)+4^'(2F),^||h 



F 

<f>'(F) + </>'{2F) 
F 



,'l|2 



-e niax{2 7 r|< ? !)"(2F)| +2eC ,4e(7r 2 +C )}||it'||| 
/or aH u GLi, where is the Lipschitz constant defined in (|5.4p . in addition, as e — > 0, 



inf / 2£ 7y'"'"U ^L-- (f) + 4,-- (2f ),gl 1 
^ew\{o} ||u'||f 2 [ 



F J 



0(e). 



Proof. Using expression (jlO.lj) for the second variation of £Q NL , we have 

A' 



^Q NL (y F )[ U , M ] = £ ^ 



(<l>"(F)+4<t>"(2F))\\P t u' t \ 



/II2 , Wll.r nUlU 



F 



-llte-P/KH' 



+ 2uJ • (</>"(2F)(P + F,_! - 2P,)K 



^(2F) 
2F 



A 

E 

£=1 
A 

E 

^=1 

/0'(2F 2 ) 0'(2F)\ 



|| (I 2 P e )(u' e+1 + u' e )\\ 2 - £ y'(2F)||P,< +1 | 



'(2F 2 )-^(2F))||P,K +1 +^)| 



V 2F 2 

N 



2F / 



||(/ 2 -P)K +1 + ^)| 12 



- E 

£=A+2 



(0"(F) + 4 0"(2F))||P^| 



/||2 , 0'(F)+20'(2F) M . T . /||2 



P 



6"(F)+4^"(2F))||Pui| 



0'(F)+2 0'(2F) 



F 



lite I 



'(F)+40"(2F))||P^ +1 f + 



^(F) + 2^(2F) 
F 



IKIa-PKjr+ill 



2e«i • U"(2F)(P 1 - P0 ui + 2e< +1 • U"(2F)(P K+1 - P A - +1 ) < +1 



Due to the assumptions 4>'(2F) > and <fi"(2F) < 0, we can drop the second sum and estimate the remaining 
terms to get 



8 2 £^ L (y F )[u,u} > min jV'(F) + 4 0"(2F), ^W|| 



|(2l) 



+ min<U"(F) + 4 0"(2F), 



II-'" 2 



+ min < 4>"{F) + 4 (j)"(2F), 



4/(F) + 2f{2F) 
F 

<j/(F) + tf/{2F)\ n ,, 



11(C) 



1^(3) 



4e 



V|| U '||| (2[) -4e 2 C (|K||| (a) + |K||| p) ) -2e7r|0"(2F)||| U '||| p) , 
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and the first part of the theorem follows. 

To obtain the sharp result for the infimum, consider again the displacements u G U and it £ IA as in the 
proof of Theorems 15. II and 15.21 Using the expansion displacement u G U, we obtain 

S 2 £^ L (y F )[u,u] = {f{F)+A<t>"{2F) + 0{e))\\u'\\%, 

where the 0(e) term arises from the interfacial terms. We note that even if the support of u was restricted 
to a smaller set, a 0(e) term would still arise due to the terms from the boundary of the support. 

On the other hand, if the support of the zig-zag displacement u G U is again chosen so that there is no 
contribution to <5 2 £^ NL from the continuum or interfacial regions, we can obtain 



where again the 0(e) term arises from the term on the boundary of the support of u. The result now 
follows. □ 

10.3. Modeling error for the quasi-nonlocal approximation. We next estimate the modeling error of 
the quasi-nonlocal approximation. Let / G 14 denote an external force applied to a deformation y F € y and 
let w a and uQ NL G U solve the linearized equations 

£ a (y F )[v}+6 2 £ a (y F )[u a ,v] = (f,v) for all v G U, 

6£ QNL (y F )[v] + 6 2 £ QNL (y F )[u QNL ,v] = (f,v) for all v G U. 

The modeling error of the quasi-nonlocal approximation, r, is again given via the duality relationship 

(t,v) :={d£^ h (y F )[v]~d£ a (y F )[v]} + {S 2 £^ h (y F )[u a ,v]~S 2 £ a (y F )[u a ,v]} for aU v G U. (10.2) 

Yet again, the solution u^ NL will play no role in the analysis of the modeling error, so to simplify the 
notation in the rest of this section we will sometimes suppress the superscript and simply write u instead of 
u a . However, in the statements of the theorems the proper notation will be used. 
We compute that we have for all v G U 

N 



5£^{y F )[v]-8£*{y F )[v]=e £ 



e=K+2 
+ e 



b'(2F 1 ) <j>'(2F 2 ) 



2F 1 
'(2F 1 ) 



2F t 
b'(2F 1 ) 



2F 2 

b'(2F 2 ) 

2F 2 
b'{2F 2 ) 



(All' + F 2 1l"' ) r -2 ^'( 2F l) III 



v't 



(2y' FA -ey' FA )+e^i^y' FA 



2F 1 

( 2 y'F,K+i + £ yF,K+ 2 ) -e 



2F, 



2F X 2F 2 
and thus observe that 

6£ ClNL (y F ) — 5£ a (y F ) =0 if y F is a linear chain since F\ = F 2 and y F = y F = 

but that 

S£^ NL (y F ) — 6£ a (y F ) ^0 if y F is a circular chain. 
After manipulations analogous to those in Section we can see that 

N 



Uf,k+ 2 



(10.3) 



J K+l 



S 2 £^ L (y F )[u,v] - 6 2 £ a (y F )[u,v] = -e ]T 

e=K+ 2 



A 



A,- A, 



£U e+1 



-u'l 



2 [A. 



■v't 



2A e j u' t 

(10.4) 



+ e 2 u" ■ A v[ - £ 2 n!' K+2 ■ A K+1 v' K+1 

+ 2eu[- (Ax - A ) v'i +2eu' K+1 ■ (a k +i - Afir+i) v' K+ i, 

where the matrices At and At are defined in (|7.3[) . We now, analogously to Theorems 17.11 and 17.21 obtain 
the following results for 1-D constrained (Theorem 110. 4p and 1-D unconstrained (Theorem 110. 5p chains. 
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Theorem 10.4. Let yp G y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe, let f G U, and let u a satisfy S£ a (yF)[v] + 5 2 £ a {yp)[u a , v] = (/, v) for all v G U. The modeling error 
of the quasi-nonlocal approximation, t, then satisfies the inequality 

1/2 



|r|U< £ |0"(2F)| 



l(u )ilr + e||(«°) 



A+2l 



N 



1/2- 



e E IK"°)/+il 



e=K+2 



Proof. Since in this case Pi = Pg = P for all £, Fx = F 2 = F, and (I 2 — P)w = for all w GU, the expression 
(|10.4[1 for the modeling error reduces to (recall that S£ (:iNh (yp) — S£ a (yp) = for linear chains) 

(r, v) = -e J2 e 2 4>"{2F){u*)r +1 ■ v' £ + e 2 cj>"{2F) (V)'/ ■ v[ - (u^ +2 ■ v' K+1 ) 



N 



1/2 



<e 2 W'{2F)\\e J2 II Will 

V t=K+2 



WWe^e.) 



+ e 



\4>"{2F)\ U||K)ri| 2 + £ ||( u % +2 || 2 \W\\ 



1/2 



1*1(3), 



and the result follows. 



□ 



Theorem 10.5. Le< yF £ y denote the 1-D configuration of atoms with nearest neighbor interatomic spacing 
Fe, let f £14, and let u a satisfy S£ a (yp)[v] + 5 2 £ a {yp)[u a , v] = (/, v) for all v G £/. TTie modeling error 
of the quasi-nonlocal approximation, t, then satisfies the inequality 

\ 1/2" 



|r||* < e max |0"(2F)| 



b'(2F) 



2F 



e\\(uXr+e\\(u a yk + 2\\ 



1/2 



e E IK«°)i+il 



=A"+2 



(10.5) 



Proof. In this case, Pf = Pi = P for all £ and Fi = i<2 = F, but we do not, in general, have (I 2 — P)w = 

for all w eU. Defining A = (j)"(2F)P+ ^—^-(I 2 — P), the expression (|10.4|) for the modeling error reduces 

2F 

to (recall that S£ ^ nl (2/_f) — S£ a (yp) = for linear chains) 

(r,v)=-e J2 e 2 (A(uJ^x)-4 + e 2 ((uJ{-Av'x-(u^ +2 -Av' K+ x) 

b'{2F) 



\4>"{2F)\, 
e max -I \4>"{2F)\ 



2F 
4>'{2F) 



( N 

e E 

V. l=K+2 



1/2 



'" II 2 



WW* 



2F 



£||(^|| 2 + £ ||(^' 



A' + 2 



1/2 



(3)' 



and the result follows. 



□ 



We now, analogously to Theorem 17. 3[ estimate the modeling error for the quasi- nonlocal approximation 
of a circular chain. In the proof, we will now need to estimate S£ QNh (y F ) - 5£ & {y F ) # 0. 

Theorem 10.6. Let yp G y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe, let f G U, and let u a G U satisfy S£ a (yp)[v] + 5 2 £ a {yF)[u a , v] = (/, v) for all v G IA. 
If (j> G C 3 (0,oo), then the modeling error of the quasi-nonlocal approximation, r, satisfies the inequality 

\ I/ 2 / N \ 1/2 / m \ V 2 1 



||r||, < (C l£ 2 + C 2 e 4 ) 



E 11 « 



?n 2 



«=A+2 



£=A+2 



e E HW+if 
s> e=K+2 / 
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+ (C 3 e + 6C e 2 + C e 3 ) 



: ) V2 +( 



1/2 



+ C K e 2 + W, 

where C\ and C'2 are defined in (|T. 10[) . C$ is the Lipschitz constant defined in (|5 .4[) . C K is defined in (|7.8[) . 
and 



C 3 = max \\</>"(2F)\ 



b'(2F) 



2F 



,2tt 



<^"(2F) 



6'(2F) 



2F 



C 3 = v/2 (2 e C (l + tt£)F + n\<j>'(2F)\). 



(10.6) 



Proo/. Using (|T03j) and the fact that \\y' e \\ = F, \\y' e '\\ = 2Fe^ 1 sin(Tre) < 2Ptt, and || = 4Fe~ 2 sin 2 (7re) < 
4Fir 2 for all ^, we first estimate similarly as in (|7.7[) 

|5^ NL (^)H -<55 a (^)M| < £ 2 C K ||t.'||^ (c) + £ 3 / 2 c a || v '||^ (3) . 

Using estimates similar to those in Theorem 17. 3i we then have 
(t,v) <e 2 fmaxl|^"(2F)|,'' 



4>'{2F) 



e 2 2ir 



<£"(2F) - - 
"(2F) - 



2F 
/(2F) 



1/2 



2F 

0'(2F) 



e max |</>"(2F)| 



2F 
4>'{2F) 



e 2tt 



>{2F)- 



2F 
>(2F) 



2F 



' ' \ t=K+2 

f4^cv 2 ) L J2 IIK)"II 2 ) WAm 

' \ e=K+2 I 

+ i2c ) L J2 II 

' \ l=K+2 

j+Q,/) (ellKy/f + slKn^lpy'V'll^p) 

+ 6c 0£ ) r e ||(u»)iii 2 + eiKu-y^+iii 2 ) 1/a ||t/||^ P j 



l«'IU§(£) 
1/2 

I 

1/2 



+ e 2 C K ||^|U i(£) +e 3/2 C 3 ||^|| li(a) , 
and the result follows. 



□ 



Remark 10.1. Notice that in all three scenarios, under the assumption that the number of interfaces between 
the atomistic and continuum regions remains constant, the modeling error is of order 0(e 3 / 2 ) as e 0. 



10.4. Error analysis for the quasi-nonlocal method. As in Section [FJ we easily see that 

(r, U a - U Q NL ) = S 2 £^ h (y F )[u a - u^ h ,u a - u^ h ], 



(10.7) 



which can be combined with the stability results and the modeling errors from the previous sections to obtain 
estimates of the errors of the quasi-nonlocal method. 

Let us first define the stability constants (see Theorems 110.11 110.21 and 110.3)) 



73 = <f>"(F) + 4 4>"{2F) and 74 = min <^ <j>"{F) + 4 (j)"{2F), 



F 

Combining f| 10 . 7[) with Thcorcms ll0.ll and ll0.41 we then obtain the following theorem for the 1-D constrained 
chain. 

Theorem 10.7. Let y F € y denote the 1-D configuration of atoms with nearest neighbor interatomic 
spacing Fe. Given f eU, let u a , u QNL £ U satisfy 5£ a (y F )[v]+5 2 £ a (y F )[u a ,v} = (f,v) and 6£ QNL {y F )[v} + 
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8 2 £ QNL {y F )[u QNL , v] = (f, v), respectively, for all veil. If <f>"(2F) < and 73 > 0, then 



73 



(e\\(uX\\ 2 + e\\(u^ +2 \\A 1/2 +s(s ]T HK^ill 2 ) 
^ ' V e=K+2 / 



1/2 



Combining (|10. 7[) with Theorems 110.21 and 110.51 we have the following theorem for the 1-D unconstrained 
chain. 

Theorem 10.8. Let y F € y denote the 1-D configuration of atoms with nearest neighbor interatomic 
spacing Fe. Given f eU, letu a , U Q NL eU satisfy 6£ a (y F )[v]+5 2 £ a (y F )[u a ,v] = (f,v) and 5£ QNL {y F )[v] + 
5 2 £ QNL (y F )[u QNL ,v] = (f,v), respectively, for all veU. If<f>'(2F) > 0, <f>"(2F) < 0, and 74 > 0, then 



\\(u a -u QNL )'\\ e 2 < Ce 



e u" 



+ 4(u a )'k + 2\\ 



1/2 



e \ e 



N \ 

E iiw+if 

=K+2 / 



1/2 



where 



C = 7 4 " 1 niax||0"(2F)|, 



b'(2F) 



2F 



Finally, note that the stability "constant" in Theorem 110.31 depends on e. To simplify the notation, we 
define 

le = 74 - £ max {2^|0" (2F) I + 2 e C 0) 4e(^ 2 + C )}, 

where is the Lipschitz constant defined in (|5.4[) . Combining (|10. 7[) with Theorems ll0.3l and ll0.6i we then 
have the following theorem for the uniform circular chain. 

Theorem 10.9. Let y F G y denote the uniform circular configuration of N atoms with nearest neighbor 
interatomic spacing Fe. Given f € U , let u a ,u^ NL G U satisfy 8£ a (y F )[v] + S 2 £ a (y F )[u a ,v] = (f,v) and 
S£ QNL (y F )[v] + S 2 £Q NL (y F )[uQ NL ,v] = (f,v), respectively, for all v e U. If <f> € C 3 (0,oo), (f>'(2F) > 0, 
4>" (2F) < 0, and i/74 > and e is small enough so that "f E > 0, then 



1 - 1 (C 1 e 2 + C 2 e 4 ) 



N 



e 



E ii(«°)*u a 



v £=K+2 



E IIM"!! 2 

1/2 



e E HWil 



fcA+2 



f=A+2 



eiiM'iir + siiMA+iir +uii(« a )iir+e|i™ +2 i 



1/2 



+ l7 \C K e 2 + C^/ 2 ), 
where C\ and C2 are defined in (|7.10[) . C3 and C3 in (|10.6p . in 



, and C K in (|7.8p . 



Remark 10.2. Stability estimates for quasi-nonlocal approximations of£ a ' b (y) that include the bond-angle 
energy, 

£ QNL - b (y), can be obtained by including the increased stability (|6.5[) . This increased stability for the 
atomistic and Cauchy-Born model was analyzed in Section® The modeling error analysis for £^ NL ' b {y) 
is identical to the modeling error analysis for £^ NL (y) since the bond-angle energy is not approximated. 
Corresponding error estimates for £Q NL ' b (y) can then be obtained from these stability and modeling error 
estimates and are summarized in the conclusion section below. 



11. Conclusion 

We summarize the results in this paper in the following tables for the case 4>'(2F) > and 4>"(2F) < 0. 
We first give results for the stability of the fully atomistic model. These stability bounds (and those that 
follow for the Cauchy-Born and quasi-nonlocal approximations) neglect 0(e) terms. We have modeled the 
strength of the bond- angle energy (|6 . 1 1) by a > 0. 
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Atomistic model 


Stability 


ID (constrained) 


cf>"{F) +4 0"(2F) > 


ID (unconstrained) 


min | (j)"(F) + 4 0"(2F), ^ + ^ } > 


circle 





We next give stabilty and error estimates for the Cauchy-Born approximation. The second-order error 
estimates 0(e 2 ) for the Cauchy-Born approximation require that ||(w a )'"||^2 be bounded uniformly in e, 
which is not the case for the approximation of atomistic configurations with defects. 



Cauchy-Born model 


Stability 


Error 


ID (constrained) 


0"(F) + 4(/)"(2F) > 


0(e 2 ) 


ID (unconstrained) 




0(s 2 ) 


circle 




0{e 2 ) 



The lack of accuracy of the Cauchy-Born approximation for problems with defects is the motivation for 
the development of atomistic-to-continuum methods such as the quasi- nonlocal method which attain 0(e 3 ^ 2 ) 
accuracy for problems with defects [3[T4j[24]. We summarize below our stability and error estimates for the 
quasi- nonlocal approximation. 



QNL model 


Stability 


Error 


ID (constrained) 


4>"{F) +4 0"(2F) > 


0(e 3 / 2 ) 


ID (unconstrained) 


min { r(F) + 4 <j>"{2F), ^p- + ^ J > 


0( e 3/2) 


circle 


f ,,, . ,„. . (J)'(F) 2acosB E ) 
min U"(F) + 4 0"(2F), Z±J. + j > 


0{e 3 ' 2 ) 
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